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ABSTRACT 

A GRB afterglow has been commonly thought to be due to continuous deceleration 
of a postburst fireball. Many analytical models have made simplifications for deceler- 
ation dynamics of the fireball and its radiation property, although they are successful 
at explaining the overall features of the observed afterglows. We here propose a model 
for a GRB afterglow in which the evolution of a postburst fireball is in an intermediate 
case between the adiabatic and highly radiative expansion. In our model, the afterglow 
is both due to the contribution of the adiabatic electrons behind the external blastwave 
of the fireball and due to the contribution of the radiative electrons. In addition, this 
model can describe evolution of the fireball from the extremely relativistic phase to the 
non-relativistic phase. Our calculations show that the fireball will go to the adiabatic 
expansion phase after about a day if the accelerated electrons are assumed to occupy the 
total internal energy. In all cases considered, the fireball will go to the mildly relativistic 
phase about 10^ seconds later, and to the non-relativistic phase after several days. These 
results imply that the relativistic adiabatic model cannot describe the deceleration dy- 
namics of the several-days-later fireball. The comparison of the calculated light curves 
with the observed results at late times may imply the presence of impulsive events or 
energy injection with much longer durations. 

Subject headings: gamma-ray : bursts — radiation mechanisms : nonthermal 
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1. INTRODUCTION 



Even though the energy source for gamma-ray bursts (GRBs) has remained unknown, 
the popular theoretical explanation for their radiative properties has been commonly 
thought to be the fireball+shock wave model, in which a GRB results from the dissipation 
of the kinetic energy of a relativistically expanding fireball. This dissipation can be either 
(most likely) due to internal shocks produced during the collision between the shells with 
different Lorentz factors in the fireball (Rees & Meszaros 1994; Paczyhski & Xu 1994; Sari 
& Piran 1997), or due to external shocks (a forward blastwave and a reverse shock) formed 
by the fireball colliding with the surrounding medium (Rees & Meszaros 1992; Meszaros 
& Rees 1993; Katz 1994; Sari, Narayan, & Piran 1996). After the main GRB, the fireball 
will continuously decelerate due to more and more swept-up medium matter and therefore 
will produce delayed emission at longer wavelengths, an afterglow, as predicted in advance 
of the observations (Paczyhski & Rhoads 1993; Katz 1994; Meszaros & Rees 1997; Vietri 
1997a). 

Afterglows from GRBs have been observed from a number of objects at X-ray, optical, 
and in one case also at radio wavelengths. The published analytical models are successful 
at explaining the major features of the light curves (Meszaros & Rees 1997; Wijers, 
Rees, & Meszaros 1997; Reichart 1997; Waxman 1997a,b; Tavani 1997; Vietri 1997a,b; 
Katz & Piran 1997; Dai & Lu 1998a). However, such models cannot provide a detailed 
description for the evolution of a postburst fireball and the light curve of an afterglow 
because they have made simplifications in three aspects. First, all of these models have 
assumed that the postburst fireball is extremely relativistic. Second, Wijers et al. (1997), 
Waxman (1997a,b), and Reichart (1997) considered the adiabatic expansion of the fireball 
in uniform interstellar medium. This is a reasonable assumption if the timescale for cooling 
of the accelerated electrons behind the blastwave is much longer than the expansion 
timescale of the fireball or the electrons carry a small fraction of the internal energy. This 
simple model has given a scaling relation between the fireball's Lorentz factor (7) and 
the blastwave's radius (R): 7 cx Dai & Lu (1998a) further discussed the effect 

of radiative corrections and nonuniformity of the medium on GRB afterglows. On the 
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other hand, Vietri (1997a,b) postulated that the postburst fireballs are highly radiative. 
This requires that the accelerated electrons behind the blastwave occupy all the internal 
energy, and that they always cool much more rapidly than the fireball expands. Such a 
model has given another scaling relation: 7 oc R"^. Third, it is usually assumed that the 
distribution of the electrons behind a shock is a power law. Since electrons with different 
Lorentz factors should have different efficiencies for synchrotron radiation in the same 
magnetic field, the spectrum radiated from higher-energy electrons for this distribution 
is steeper than that from lower-energy electrons (Sari, Piran, & Narayan 1998). The 
adiabatic expansion model (Wijers et al. 1997; Waxman 1997a,b; Reichart 1997; Dai 
& Lu 1998a) has assumed all of the accelerated electrons behind the blastwave to be 
adiabatic, while in the radiative expansion model (Vietri 1997b) these electrons has been 
thought to be highly radiative. 

In fact, the fireball first expands relativistically, and will eventually go into the non- 
ralativistic phase (the Scdov phase) after a long time. Furthermore, the actual expansion 
of the fireball is likely in an intermediate case between the adiabatic and highly radiative 
expansion. Finally, an afterglow may be contributed by both the adiabatic electrons 
and radiative electrons behind the blastwave. In this paper we would propose a model 
which addresses the above issues. Huang et al. (1998) first studied numerically the 
evolution of an adiabatic fireball from the ultrarelativistic expansion phase to mildly 
relativistic expansion phase. The present work is in fact a significant development of the 
study of Huang et al. (1998) through considering the above three issues. This paper 
is organized as follows: in section 2 we calculate numerically the dynamical evolution 
of a postburst fireball from the ultrarelativistic expansion phase to Sedov phase. In 
section 3 we formulate synchrotron radiation from the accelerated electrons behind the 
blastwave and compare our results with observations, and in the final section wc give a 
brief discussion. 

2. HYDRODYNAMICS OF POSTBURST FIREBALLS 
We assume that a fireball with an amount of energy E comparable to that observed 



3 



in gamma rays, E ~ 10^^-10^^ ergs, and with the mass of the contaminating baryons, Mq, 
is produced, and after an initial acceleration its Lorentz factor is 77 = E/{Mq(?). Subse- 
quently, at the radius Rq, the expansion of the fireball starts to be significantly influenced 
by the swept-up medium and two external shocks may form: a forward blastwave and a 
reverse shock (Rees & Meszaros 1992). As usual, Rq is supposed to be 

where £^51 = £^/10^^ergs, 77300 = 7^/300, and tio is the electron number density of the 
medium in units of lcm~^ (tio = 7i/lcm~^). Following the main GRB event, which may 
be produced by nonthermal processes such as synchrotron or possibly inverse Compton 
emission, the blastwave continues to sweep up the medium. 

According to Blandford & McKee (1976), the electron number density {n') and energy 
density (e') of the shocked medium in the frame comoving with the fireball and the Lorentz 
factor of the blastwave (F) can be written as 

= (2) 

7-1 

e' = II±i(7-l)7i7n,c2, (3) 
7-1 

r [ (7 + l)[7(7-l) + lF ]'^' m 
1 7(2-7)(7-l) + 2 j ' ^ ^ 

where 7 is the adiabatic index of the shocked medium, which is generally between 4/3 
and 5/3. One expects that equations (2)-(4) are appropriate for both relativistic and 
non-relativistic blastwaves. From the definition of 7 (Blandford & McKee 1976), we have 

derived an expression: 7 ^ (47 + l)/(37). It can be seen from this approximation that 
7 ~ 4/3 for an extremely relativistic blastwave and 7 5/3 for the Sedov shock. 

We further assume that the magnetic energy density in the comoving frame is a 
fraction of the total thermal energy density, viz., B' — ^^(STre')^/^, and that the 
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accelerated electrons behind the blastwave carry a fraction of the energy. This implies 
that the minimum Lorentz factor of the random motion of electrons in the comoving frame 
is 7mm = ^e{l — i)'>Tip/'me + 1. We here consider only synchrotron emission from these 
electrons, and neglect the contribution of inverse Compton emission because the latter 
emission is not of importance particularly at late times of the evolution (Waxman 1997a; 
Dai & Lu 1998a). The energy of a typical accelerated electron behind the blastwave 
is lost both through synchrotron radiation and through expansion of the fireball, and 
thus the radiative efficiency of this electron is given by t'~l^/{t'~l^ + f'J-), where t'^y^ is 
the synchrotron cooling time, t^^^ = 67rmec/((Tr-B'^7^j„), and t'^^ is the comoving- frame 
expansion time, t'^^ = R/i'yc) (Dai & Lu 1998a). Here R is the radius of the blastwave. 
Since all of the accelerated electrons behind the blastwave carry only a fraction of the 
internal energy, the radiative efficiency of the fireball can be given by 



For the adiabatic expansion, ^ 1 or ^ t'^^, so / ~ 0; but for the highly radiative 
expansion, ~ 1 and t'^y^ -C lead to / ~ 1. One expects that in the intermediate 
case the radiative efficiency of the fireball (/) is between and 1. 

In the absence of radiation and expansion losses, according to equations (2) and (3), 
the kinetic energy per baryon of the shocked medium in the comoving frame is (7 — l)mpC^. 
But, in the presence of these losses, such a kinetic energy becomes (7 — 1)(1 — f)mpC^. 
Thus, the total kinetic energy of the shocked medium in the burster's rest frame is {7[(7 — 
1)(1 — /)-|-l] — l}Mc^, where M is the mass of the swept-up medium, M — {A7r/3)R^nmp. 
We therefore obtain the total kinetic energy of the fireball 



where the second term is the kinetic energy of the contaminating baryons. If this term 
is neglected and / w 0, then Ek ~ j^Mc^ for extremely relativistic expansion. This 

expression is just the starting point of many works (e.g., Mcszaros & Rccs 1997; Wijers 
et al. 1997; Waxman 1997a,b; Dai & Lu 1998a). lif^l, then w (7 - 1)(M + Mo)c2, 




(5) 



Ek = {7[(7 - 1)(1 - /) + 1] - IjMc^ + (7 - l)Moc2 , 



(6) 
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which turns out to be the case discussed by Vietri (1997a,b) and Katz & Piran (1997). 
It should be pointed out that at late times the kinetic energy of the fireball is ~ 
(7^ — l)Mc^, which is a factor of 1.4 larger than the Sedov result (Blandford & McKee 
1976). Hence, equation (6) can also describe well the non-relativistic evolution of the 
fireball. Due to synchrotron radiation, the rate for the kinetic energy loss is given by 
(Blandford & McKee 1976) 

— ^ = -4rf(/3c)7(7 - l)nm,cV , (7) 

where tb is the burster-rest-frame time and /3 = (1 — 7~^)^^^. 

In order to study dynamical evolution of the fireball, we should add two differential 
equations (Huang et al. 1998): 

^ = (i-r-)V^c, (8) 

^ = 7(7 + \/7^), (9) 

where t is the observer-frame time. 

Equations (3)-(9) present a perfect description for the dynamical evolution of the 
postburst fireball. In order to solve these equations, we must determine the initial condi- 
tions. Assuming that 70 and E^o are the initial values of the Lorentz factor and kinetic 
energy of the fireball respectively, we require 

^feo = I = {7o[(7o - 1)(1 - /o) + 1] - l}(MocVr^) + (70 - l)Moc^ , (10) 

where /o is the initial value of the radiative efficiency. Let's define an index a through 
the following expression: 

din 7 

Q; = — n — ■ (11) 

We take E = 10^^ ergs, n = 1 cm~^, and Mq = 10~^ Mq and 2x lO^^M©, and our numerical 
results are shown in Figures 1-7. Figure 1 illustrates evolution of the fireball's kinetic 
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energy, and Figures 2 and 3 give evolution of j{t) and R{t). Figure 2 shows that the 
fireballs will go to the mildly relativistic phase about 10^ seconds later, and to the non- 
relativistic phase after several days, implying that the relativistic adiabatic model cannot 
describe the deceleration dynamics of a several-days-later fireball. We have compared 
our numerical result with the analytical solution for the Sedov phase, and found that 
they are in good agreement at the non-relativistic phase. This comparison is shown in 
Figure 4. In Figures 5 and 6 the time-dependence of / and a is plotted. It can be seen 
from Figure 5 that in the case of = 1 (the solid, dotted and dashed lines) / is first 
kept to be a constant (~ 1), subsequently declines quickly, and finally tends to zero after 
10^ seconds, showing that the fireball first expands radiatively, soon later goes into the 
intermediate expansion phase, and finally becomes a non-relativistic shock. In this case, 
a first increases up to a peak near 3 due to the influence of the contaminating baryons, 
and then decreases to 1.5 at about 1 day. After this, a decreases to zero because the 
Lorentz factor becomes one during the non-relativistic phase. Figures 7 gives evolution 
of the mass of the swept-up medium. 



3. X-RAY AND OPTICAL RADIATION 



In the absence of radiation loss, the distribution of the accelerated electrons behind 
the blastwave is usually assumed to be a power-law function of electron energy: 

^ OC 7~^, for 7^j„ < 7e < 7max , (12) 

where 7^,100; is the maximum Lorentz factor, -f^nax = 10^{B' /IG)-^/^ (Dai & Lu 1998a), 
and p is the index between 2 and 3. However, radiation loss may modify such a simple 
distribution. In a magnetic field, electrons with different Lorentz factors have different 
efficiencies for synchrotron radiation. As defined by Sari et al. (1998), the critical electron 
Lorentz factor, 7c, above which synchrotron radiation is significant, is written as 

^ ^ 3me 1 .^g. 
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Electrons with Lorentz factors below 7c are referred to as adiabatic ones, and electrons 
above 7c as radiative ones. In the presence of steady injection of electrons accelerated by 
the shock, the distribution of radiative electrons becomes another power-law function with 
an index of p + 1 (Rybicki & Lightman 1979), but the distribution of adiabatic electrons 
is unchanged. Thus, the actual distribution can be given in three cases: (i) For 7c < 7mm, 

^^C,^;^P+^\ C,= _^ Ntot, (7mm < 7e < 7ma.) , (14) 

where N^ot is the total electron number of the shocked medium {Ntot = M/rUp). (ii) For 



d^e ^ f C27e"^ > 7min < 7e < 7c 

d^e I Cs'y-^P+^^ , 7c < 7e < Imax 

where 



(15) 



C2 = C3/7c, (16) 

Ntof (17) 



C: 



I mm ic _j_ I max 



lc{p - 1) P 

(iii) If 7c > 7maa;, then 

dNe _ P — 1 

-J-^ = Ci%P , C4 = r^Nt^t , (7mm < 7e < 7max) • (18) 

"7e 7mm ~ "Ymax 

After having the modified electron distribution functions, we can calculate the radi- 
ation flux. The power for synchrotron radiation from all the accelerated electrons of the 
shocked medium in the comoving frame is given by 



where 



r+00 

F{x) = X / K5/3{t)dt, (20) 

Jx 



8 



and 



lleB' 



(21) 



27rmeC 



with K^i^it) being the Bessel function. The observer-frame flux density should be 



5. = 7(l + /3fj 



7(1 + 13)) AnD 



i2 ' 



(22) 



where D is the source distance to the observer. In writing equation (22), we have assumed 
that the emitting equal-time surface of the source is an ellipsoid. However, the actual 
fireball always decelerates due to more and more swept-up medium. It has been found 
that due to the deceleration the emitting surface becomes a distorted ellipsoid (egg-like 
shape) (Panaitescu & Meszaros 1998; Sari 1998), which slightly influences the observer- 
frame flux density. Huang et al. (1998) considered this effect by introducing a factor near 
five. Here we would neglect this effect. The flux observed by an X-ray detector is an 
integral of S,y: 



where Uu and ui are the upper and lower frequency limits of the detector. 

In the previous section, the dynamical evolution of a postburst fireball has been 
evaluated numerically. Now we continue to calculate the afterglows at X-ray and optical 
wavelengths. While some of the parameters are fixedly taken (E — 10^^ ergs, n — 1 
cm~^, p = 2.1), we change the others such as Mq, ^'^ and in reasonable ranges so as 
to investigate their influence. The results are plotted and compared with observations 
in Figures 8 and 9. In these figures, the solid lines correspond to Mq — 2 x 10~^ Mq, 
^1 = 0.01, = 1; the dotted lines correspond to Mq = 1 x 10~^ Mq, ^| = 0.01, = 1; 
the dashed hues correspond to Mq = 2 x 10~^ Mq, = 0.1, — 1; and the dash-dotted 
hues correspond to Mq = 2 x 10~^ Mq, ^| = 0.01, = 0.5. Figure 8a is the R band 
afterglow from GRB 970228 and Figure 8b illustrates the optical afterglow from GRB 
970508. Although we have assumed different distances in Figures 8a and 8b, it is not 




(23) 
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necessary to imply that GRB 970228 lies farther from us than GRB 970508 does, since 
other intrinsic parameters such as n may be different. X-ray afterglows or initial X- 
ray bursts were detected for nine GRBs. Table 1 hsts all the flux data available in the 
hterature. In Figure 9, X-ray afterglows from these GRBs are plotted together. Please 

note that since different detectors work in different bands, here we have converted the 
flux data into 0.1 — 10 keV band linearly, errors to a factor of two are thus possible. 

We see that the present model generally fits the observations and small changes of 
parameters do not alter the overall properties of afterglows. However, a problem appears 
at later times {t > 10^ s) after the fireball becomes non-relativistic. The model predicts a 
sharper decline while the observed fiuxes are obviously higher. In fact, the optical after- 
glow from GRB 970228 was observed to follow approximately a power-law decay for at 
least ~ 190 days (Pruchter, Bergeron, & Plan 1997) and the afterglow fi:om GRB 970508 
decayed even more slowly after about 80 days, implying the presence of a constant compo- 
nent (Pedersen et al. 1998; Galama et al. 1998a). The overall power-law decay lasting for 
several months is usually considered as strong evidence for the fireball+blastwave model. 
Here we would like to stress that the blastwave in a simple fireball model will cease to be 
relativistic after ~ 10^ s. At least, it is problematic to assume a blastwave with 7 > 5 
when t > 10^ s (also see Huang et al. 1998). One should be cautious in applying the sim- 
ple scahng laws such as 7 oc R oc t^^^, and Si, oc in the adiabatic expansion 
model at such late times. So the long-term optical afterglows have really raised a problem 
to the popular fireball-|-blastwave model. For GRB 970508, the optical fiux peaked about 
two days later, possibly associated with an X-ray outburst (Piro et al. 1998). This feature 
could not be explained by a simple fireball-|-blastwave model, too. Our theoretical light 
curve peaks several hours later at optical wavelength and several tens of seconds later in 
X-rays. 

4. DISCUSSION 

In this paper we have tried to propose a model for a GRB afterglow in which the 
evolution of a postburst is in an intermediate case between the adiabatic and highly 
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radiative expansion. In this model the afterglow is due to the contributions both of the 
adiabatic electrons and of the radiative electrons behind the blastwave. In addition, our 
model is valid both for the extremely relativistic phase and for the Sedov phase. Our 
calculations show that the postburst fireball will go to the adiabatic expansion phase 
after about a day if the accelerated electrons arc assumed to occupy the total internal 
energy (viz., = !)• In all cases considered, the fireball will go to the mildly relativistic 
phase about 10^ seconds later, and to the non-relativistic phase after several days. These 
results imply that the relativistic adiabatic model (7 oc t~^^^ and R oc t^^^) isn't suitable 
to describe detailedly the deceleration dynamics of the several-days-later fireball. What 
we would like to emphasize is that one should be cautious in applying the simple scaling 
laws such as 7 oc t~^^^, R oc t^^"^, and Si, oc t3(i-p)/4 ^j^g adiabatic expansion model at 
t >10^ seconds. 

In comparing our model with observations, we find that for GRB 970228 the initial 
fireball can produce the observed optical afterglow, but the flux of the flrcball after several 
days is below the observed data, since the observed afterglow follows approximately a 
power-law decay for at least ~ 190 days. The optical afterglow from GRB 970508 first 
declined slowly, then rised to a peak at about two days and after this it decayed. But the 
light curve became fiatter after 80 days than the power-law decay. Comparison of these 
observational results with our calculations may imply the presence of impulsive events or 
energy injection with much longer durations. One possibility is that a postburst fireball 
contains shells with a continuous distribution of Lorentz factors (Rees & Meszaros 1998). 
As the external blastwave sweeps up ambient matter and decelerates, internal shells will 
eventually catch up with the blastwave and supply energy into it. A detailed calculation 
shows that this model can explain well the afterglow from GRB 970508 (cf. Panaitescu, 
Meszaros & Rees 1998) 

Another possibility is that a strongly magnetized millisecond pulsar can supply its 
rotational energy into a postburst fireball (Dai & Lu 1998b). Many energy-source models 
of GRBs all predicted that as an extremely relativistic fireball is produced, a strongly 
magnetized millisecond pulsar is born. It is natural to expect that magnetic dipole ra- 
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diation from the pulsar may influence evolution of the external fireball because the elec- 
tromagnetic waves are always absorbed by the shocked medium. Such an effect has been 
analytically shown to be able to provide a satisfactary explanation for the flattening be- 
havior of the hght curve of the optical afterglow from GRB 970228 (Dai & Lu 1998b). A 
recent analysis (Dai & Lu 1998c) further shows that this effect can also explain well the 
decline-rise-decline feature of the light curve of the optical afterglow from GRB 970508 if 
the index (p) of the power-law distribution of the accelerated electrons behind the blast- 
wave increases from initial p — 1 to p — 2.2 two days later, inferred from the observed 
spectrum (Galama et al. 1998a,b). 

This work was supported by the National Natural Science Foundation of China. 
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Table 1. X-ray observations of recently localized GRBs. 
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FIGURE CAPTIONS 



Figure 1. Evolution of tlie fireball's kinetic energy. The solid lines is drawn with "stan- 
dard" parameters (E = 10^^ ergs, n = 1 cm^^ Mq = 2 x 10^^ M©, ^% = 0.01, = 1-0). 
Each of the other lines is plotted with one parameter altered. The dotted line corresponds 
to Mo = 1 X 10^^ Mq, the dashed line corresponds to = 0.1, and the dash-dotted line 
corresponds to .^e = 0.5. 

Figure 2. Evolution of the fireball's Lorentz factor. Parameters are the same as in Figure 
1. 

Figure 3. Evolution of the fireball's radius. Parameters are the same as in Figure 1. 

Figure 4. Velocity of the blastwave (log V/ c) vs. the fireball's radius. Full line is our 
numerical result, plotted with the "standard" parameters as in Figure 1. Dotted fine is 
the analytic result for non-relativistic phase, = 1.15^ x AE/{25nV'^) (Lang 1980; Shu 
1992), in which we have set E — 7.2 x 10^^ ergs, since this is approximately the kinetic 
energy left in a non-relativistic fireball for our model (see Figure 1). 

Figure 5. Evolution of the radiative efficiency of the fireball. Parameters are the same 
as in Figure 1. 

Figure 6. Evolution of the index a, where a is defined as 7 oc R~°'. Parameters are the 
same as in Figure 1. 

Figure 7. Evolution of the mass of the swept-up medium. Parameters are the same as 
in Figure 1. 

Figure 8. Predicted optical afterglows in R band. Sr. We take p = 2.1, D = 3 Gpc (a) 
or D = 1.5 Gpc (b). Other parameters are the same as in Figure 1. Also plotted are the 
observed afterglows from GRB 970228 (a) (Galama ct al. 1998a; Fruchter et al. 1997; 
Garcia et al. 1997; Wijers, Rees, & Meszaros 1997), from GRB 970508 (b) (Pedersen et 
al. 1998; Garcia et al. 1997; Galama et al. 1998a). 

Figure 9. Predicted X-ray afterglows (0.1-10 keV). The fiux is in unit of ergs cm~^ s~^. 
We take D = 3 Gpc. Other parameters are the same as in Figure 8. Also plotted are the 
observed data, which have been linearly scaled to 0.1 — 10 keV fiuxes from Table 1. 
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